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Abstract 

The Nambu Jona-Lasinio model with a Polyakov loop is extended to finite isospin chemi- 
cal potential case, which is characterized by simultaneous coupling of pion condensate, chiral 
condensate and Polyakov loop. The pion condensate, chiral condensate and the Polyakov 
loop as functions of temperature and isospin chemical potential are investigated by mini- 
mizing the thermodynamic potential of the system. The resulting (T, fj,i) phase diagram 
is studied with emphasis on the critical point and Polyakov loop dynamics. The tricritical 
point for the pion superfluidity phase transition is confirmed and the phase transition for 
isospin symmetry restoration in high isospin chemical potential region perfectly coincides 
with the crossover phase transition for Polyakov loop. These results are in agreement with 
the Lattice QCD data. 



PACS number(s): 12.38.Aw; 11.30.RD; 12.38.Lg; 



1 Introduction 



Confinement and chiral symmetry breaking are the fundamental properties of QCD in the 
nonperturbative domain. In principle, the color liberation associated with deconfinement is 
a phenomenon distinguishable from the chiral symmetry restoration phase transition, but it 
is likely that the mechanism of confinement would be closely related to chiral dynamics [1]. 
It has been argued that the confined vacuum necessarily breaks the chiral symmetry and the 
energy scale of confinement is less then the energy scale of chiral symmetry breaking: Aqcd for 

E-mail: zhaozhang@pku.edu.cn 
2 E-mail: yxliu@pku.edu.cn 



1 



confinement and ~ 47r/ 7r (/ 7r is the pion decay constant) for chiral symmetry breaking [2, 3]. In 
contrast with the situation at zero temperature, the lattice data at finite temperature suggest 
that the deconfinement phase transition and chiral symmetry phase transition occur at the same 
temperature [4-9]. Despite the attempts to explain above behaviors, the underlying reason is 
still unknown. 

QCD phase diagram and thermodynamics has been the subject of intense investigation in re- 
cent years. Lattice simulations is a principle tool to explore the qualitative features of strongly 
interacting matter and make quantitative prediction of its properties. Over the years, this for- 
mulation has given us a wealth of information about the phase diagram and thermodynamics 
at finite temperature and limited chemical potential [10-18]. In response to the lattice simula- 
tion, many phenomenological models in terms of effective freedom degrees have been proposed 
to give an interpretation of the available lattice data and to make prediction in the region of 
phase diagram that can't be reached on the lattice yet. 

As an effective chiral filed theory, classical Nambu-Jona-Lasinio model (NJL) [19-21] can 
illustrate the transmutation of originally light quark or massless quark into massive quansipar- 
ticles and generate the pions as GoldStone bosons of spontaneously broken chiral symmetry at 
the same time. This type models have also been used extensively to explore color supercon- 
ductivity phase transition at high baryon density [22]. However, NJL models have a principle 
deficiency. The reduction of the local color symmetry to the global color symmetry leads to 
quarks not confined in standard NJL. To cure this problem, an improved field theoretical quasi- 
particle model, a synthesis of Polyakov loop dynamics with NJL(PNJL), has been proposed to 
interpret the Lattice QCD results and extrapolate into the regions not accessible by lattice sim- 
ulations [23-28]. For the chiral symmetry breaking, chiral condensate can be used as an order 
parameter for chiral phase transition in the chiral limit. In contrast, the deconfinement phase 
transition has a definite meaning only at an infinite quark mass, which the expectation value 
of Polyakov loop serves as an order parameter. Since the Z(3) center symmetry is explicitly 
broken in the real world with finite quark mass, no rigorous order parameter is established for 
deconfinement phase transition [29]. However, the Polyakov loop can still act as an indicator 
of a rapid crossover towards deconfinement for finite quark mass. The motivation behind the 
PNJL model is to unify the aspects of chiral symmetry breaking and confinement through intro- 
ducing both the chiral condensate and the Polyakov loop as classical, homogeneous fields which 
simultaneously couple to quarks according the symmetry and symmetry breaking patterns of 
QCD. 

The effectiveness of PNJL model has been tested in the literature by confronting the PNJL 
results with the Lattice QCD data. It has been reported that the two flavor PNJL model can 
reproduce the result that the crossovers for deconfinement phase transition and the chiral phase 
transition almost coincide[25, 28]. For finite quark chemical potential, further investigations also 
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suggest that the thermodynamics and susceptibilities obtained in PNJL model are excellently in 
agreement with the corresponding Lattice QCD data[28, 30, 32]. In addition, the PNJL model 
has been extended to include the diquark degrees of freedom and used to explore the phase 
diagrams in high baryon chemical potential[33]. Though the PNJL results have a satisfactory 
agreement with some lattice data, there still need to perform further test of the PNJL since 
for finite baryon chemica potential the available lattice simulation is limited only in the low 
region. 

It is well known that there is no sigh problem for the lattice calculation for the case with finite 
isospin chemical potential and zero baryon chemical potential [34-3 7]. Therefore, it is interesting 
to investigate the phase diagram and thermodynamics in the framework of PNJL model at finite 
isospin chemical potential and compare the results with the corresponding lattice data. In [38], 
the authors have extended the PNJL to the case with low baryon and isospin chemical potential 
without considering the pion condensate. Both the QCD effective filed theory[34] and lattice 
calculations [35-37] have shown that the pion superfluidity phase (charged pion Boson-Einstein 
condensate) may occur in high enough isospin chemical potential. This phase had also been 
explored by other chiral models of QCD such as Ladder QCD[39], NJL model[40-42], random 
matrix model[43] and global color model[45]. Therefore, it is necessary to extend the PNJL 
model to high isospin density region to investigate the influence of Polyakov loop dynamics on 
the phase structure by including the pion condensate degrees of freedom and give further test 
of the validity of the PNJL by confronting the corresponding Lattice QCD data. In addition, 
due to the limitation of lattice simulation, the talk between effective fields theory and Lattice 
QCD is needed. 

The primary purposes of this paper is to derive the formulation of the PNJL model with 
simultaneously considering the pion condensate, chiral condensate and Polyakov loop degrees 
of freedom at mean filed level and then use it to explore the phase diagram of two flavor QCD 
to compare with the corresponding lattice data. The emphases will be put on the critical point 
and the impact of Polyakov loop dynamics. For simplicity, the freedom degrees of diquark is 
not included in this paper and the baryon chemical potential is limited in low region. 

2 PNJL model at finite isospin chemical potential 

Extending the work [28] to finite baryon chemical potential and isospin chemical potential case, 
the lagrangian of two-flavor PNJL model is given by 

Cpnjl = ^ (il^ + l0 jx - mo - ^ 75 ri) ip + G [(^) 2 + (^75^) 2 ] - U {$[A}, *[A],T) , (1) 
where ip = (ip u ,ipd) T is the quark field, 

D f * = d»-iA f * and A* = S^A . (2) 
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The two-flavor current quark mass matrix is tjiq = diag(m u ,md) and we shall work in the 
isospin symmetric limit with m u = = mo. The quark chemical potential matrix fi takes the 
form 



with 



/i 



3 



— — and 



Pi 



H- m 



— Pd 



(3) 



(4) 



2 3 r 2 ' 

where is baryon chemical potential corresponding to conserved baryon charge and m is 
isospin chemical potential corresponding to conserved isospin charge. A term proportional to A 
is introduced in (1) which explicit break the t// 3 (l) symmetry (Below we call it I3 symmetry). 
A local, chirally symmetric scalar-pseudoscalar four-point interaction of the quark fields is 
introduced with an effective coupling strength G. 

In comparison with the classical NJL lagrangian, the gauge field term A fM (x) = gAa{x)^f 
is contained and an effective potential T) expressed in terms of the traced Polyakov 

loop $ = (Tr c L)/N c and its (charge) conjugate <J> = (Tr c Zj)/iV c in the PNJL lagrangian. The 
Polyakov loop L is a matrix in colour space explicitly given by 



r 

L (x) = V exp i dr A4 (x, r 

L Jo 



(5) 



with (3 = 1/T the inverse temperature and A4 = iA°. In a convenient gauge (the so-called 
Polyakov gauge), the Polyakov loop matrix can be given a diagonal representation [29]. The 
coupling between Polyakov loop and quarks is uniquely determined by the covariant derivative 
Dp in the PNJL Lagrangian (1). For simplicity, the temporal component of Euclidean gauge 
field A4 is treated as a constant in PNJL, and the Polyakov loop is reduced as 



L 



V exp 



A A dr 



exp 



iAa 



(6) 



Corresponding to above expression, the trace of the Polyakov loop, <!>, and its conjugate, 
are treated as classical field variables in PNJL. 

The temperature dependent effective potential U($>, <£, T) is used to mimic pure-gauge Lattice 
QCD data, which should have exact Z(3) center symmetry. In this paper, we will use the 
potential 14(&,&,T) proposed in [28] , which takes the form 



W($,$,T) 



with 



rp4 



b2 (T) 



2 6 v 



$3) + ^ ^ 
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A precision fit of the coefficients Oj, b\ is performed to reproduce the lattice data and To = 
270MeV is adopted in this paper. 

To describe the spontaneous chiral symmetry breaking and I3 symmetry breaking, we define 
the chiral condensate as 

= cr = a u + a d , (9) 
with a u = (uu) and ad = (dd), and the charged pion condensates 

(^7 5 T+V>) = 7T + = ^e if? , {ipi^T-ip) = vr" = ^| e ~^> ( 10 ) 

where r± = (ri ± T2)j\pl and Tj is the Pauli matrix in flavor space. In the chiral limit, 
nonzero condensate a indicates spontaneous chiral symmetry breaking, and nonzero condensate 
7r indicates spontaneous I3 symmetry breaking. The phase factor 6 describes the direction of 
the I3 symmetry breaking. For convenience, we adopt 9 = in the following and the pion 
condensate can be expressed as 

(^75TlV ; ) = TT. (11) 

This choice of the 1$ symmetry breaking direction is consistent with the form of the explicit 
isospin breaking term introduced in (1). 

Using the standard bosonization techniques, the mean filed lagrangian of PNJL takes the 
form 

£mf = 4> (i-TuD* + 70A - M - Ni-ysn) ^ - G [a 2 + vr 2 ] - U ($ [A] , $ [A] , T) , (12) 
with the gaps 

M = m - 2Ga, (13) 
N = A - 2Gvr. (14) 

The thermal dynamical potential in the mean filed level is expressed as 
n = U(<S>,$,T)+G(cj 2 +TT 2 )-^ F lndetS m l 



(15) 



where the sum is taken over Matsubara frequencies w n = (2n + l)7rT and the trace is taken over 
color, flavor and Dirac indices. The momentum dependent inverse quark propagator matrix 
including both chiral and pion condensates in flavor space takes the form 

(iw n + fi + hi + iA 4 )j - 7 • p- M -i^N 

-Z75./V (iw n + 11 - 111 + ^4)70 - 7 ■ p - M 

(16) 



Sj(iw n ,p) = 
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Using the identity TV ln(X) = In det(X) and the technique 



det 



A B 
C D 



&et(A)&et(B)&et(C)det{C- l DB~ l - A 



(17) 



the thermal dynamical potential(15) can be further expressed as 



n 



G(a 2 + vr 2 ) + W(d>, T) - 2N C J ^ [E~ + E+] 6 (A 2 - f 



-2T 



d 3 p 



3 TV C 



ln[l + L + e~^ + ^ T ] + ln[l + Le^ E '-^ T ] 



+ln[l + L + e~^ + ^ T ] + ln[l + Le~^~^/ T ] 



(18) 



with quasiparticle energy Ep = y/(E p ± m) 2 + N 2 and E p = yjp 2 + M 2 . Performing the 
remaining color trace, the integrand of last term on r.h.s. (18) is given by 



In det 
+ In det 



1 + Le -( E p-»)/ T 



+ In det 



1 + Le ~(E+-n)/T 



+ In det 



1 + L + e -(E-+^)/T 
l + L + e -{E++„)/T 



= In 

+ In 

+ In 

+ In 



1 + 3 + $ e -(Ep~n)/T^ e -(E~-n)/T + e -3(E~-n)/T 

1 + 3 + $e"( B P +M ) /T ) e"(^ +At ) /T + e - 3 (^+^) /T 
1 + 3 ($ + $ e -«-/*)/T^ e -(^ + -^)/ T + e - 3 (^ + -^)/ T " 

1 + 3 + $ e -(^ + +M)/T^j e -(S P ++M)/T + e -3(E++^)/T 



(19) 



From (18) and (19), we see that the trace of the Polyakov loop $ and its conjugate <3> can still be 
factored out despite the appearance of the off-diagonal terms in the inverse quark propagator. 
It is easy seen from (19) that the classical NJL model is restored from PNJL when <3> and $ 
approach 1. 

Minimizing the thermal dynamical potential (18), the motion equations for the mean fields 
a, 7T, <!> and $ are determined through the coupled equations 



on _ on _ dn _ on _ 
fo~ 0, d^~ 0, 5¥~ ' dl -0 ' 



(20) 



This set of equations is then solved for the fields a, it, $ and $ as functions of the temper- 
ature T, baryon chemical potentials fiB and isospin chemical potential fij. Note that when 
there exist multi roots of these coupled equations, the solution corresponding to the minimal 
thermodynamical potential is favored. 

The NJL part of the model involves three parameters: the current quark mass of u and d, the 
local four fermion coupling constant G and the three- momentum cutoff A. In this work, these 
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parameter are fixed to reproduce three physical quantities in vacuum: pion mass m n = 140MeV, 
pion decay constant f w = 93MeV and chiral condensate (uu) = (dd) = cjq/2 = — (250MeV) 3 
with 

m = 5.5MeV, G = 5.04GeV" 2 , A = 0.651GeV. (21) 

Once £1 is known, the thermodynamic functions that measure the bulk properties of matter 
can be obtained. For an infinite system, the pressure p, the entropy density s, the baryon 
number density n#, the isospin number density nj, the flavor number density n u and n^, the 
energy density e and the specific heat c are all derived from 17. 

For = 0, the trace of Polyakov loop $ is equal to its conjugate <& and both fields are real. 
Without concerning the isospin chemical potential, numerical results suggests that the crossover 
phase transitions for the chiral condensate a and for the Polyakov loop <I> almost coincide 
in PNJL and the difference between the two transition temperatures is within 5MeV[28, 30]. 
Following [28], we define T c as the average of the two transition temperature and we verified 
the result that Tj? = 227MeV [30] for hb = Hi = 0. Since the PNJL model uses a simple 
static background field representing the Polyakov loop and reduce the gluon dynamics to chiral 
point couplings between quarks, this model can not be expected to work beyond a limited 
temperature and quark chemical potential. It is assumed this model works up to an upper 
limit of temperature with T < (2 — 3)T C ° because the transverse gluon degrees of freedom will 
be important for T > 2.5T C °[31]. In this paper we only consider the case for = and we 
assume that PNJL works well in the region |^/| < 500MeV 

3 Results 

3.1 For A = and fi B = 

For A = 0, the charged pion condensate 7r is the true order parameter for the pion superfluidity 
phase transition. Previous investigations based on the QCD effective field theory[34] and lattice 
simulations [35-37] have shown that pion superfluidity phase occurs when /xj > /x c = m w /2 for 
T = and [ib = 0, where the phase transition is second order. At sufficient high temperature 
and [ii > [i c , the pion condensate evaporates and the spontaneous breaking ^3 symmetry and 
parity are restored. The lattice calculations also suggest that the pion condensate vanishes 
for sufficient high [i\ and this ^3 symmetry restoration phase transition may be first order for 
large enough /ij . Therefore, there should exist a tricritical point (TCP) in the (T,fij) phase 
diagram. We will show below that these conclusions can be verified in our calculations within 
the PNJL model. 

Fig. 1 shows the temperature dependence of the pion condensate, chiral condensate and the 
trace of Polyakov loop at different fij. In comparison with the zero (ij case, the magnitude 
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of chiral condensate is suppressed by the isospin chemical and it almost keeps as a constant 
until the pion condensate evaporates at high T' c . The dependence of Polyakov loop <E> on T is 
similar to the zero \ij case. In the low \xj region, the I3 symmetry restoration phase transition 
is second order. With increasing fii, the evaporation of the pion condensate gets more and 
more abrupt and the I3 symmetry restoration phase transition eventually becomes first order. 
Fig. 1(c) shows that <E> is also discontinuous at the first order phase transition point T' c . This 
indicates there exists a TCP in the (T, fij) phase diagram, which was first proposed in [36] by 
using lattice simulation. Fig.l also shows that the I3 symmetry restoration phase transition 
temperature T c is always less then T®. 

Although chiral condensate is suppressed by the pion condensate, it still keeps considerable 
value in the low fij region with fij > fi c . Fig. 2 shows that the crossover phase transitions 
for the chiral condensate a and for the Polyakov loop also perfectly coincide at T c with 
T' c < T c < T c °. Therefore, if we define the peak of the Polyakov loop susceptibility 1 as the 
deconfinement phase transition point, the pion superfluidity phase is always confined at low jij 
region. For high fj,i, chiral condensate is greatly reduced and it has no meaning to compare the 
chiral susceptibility with the Polyakov loop susceptibility. Numerical results indicate that the 
peak of Polyakov loop susceptibility excellently coincides with the pion condensate vanishing 
point for high enough fij, while for the not too large /U/ the difference of these two points is 
within a few MeV. This points will be illustrated more clearly in the next section by considering 
a small I3 symmetry breaking term. Therefore, we get conclusion that the pion superfluidity 
phase is always confined. 

The isospin chemical dependence of a, it and are plotted in Fig. 3. For a fixed T, 
the Polyakov loop $ is always a monotonic increasing function of m and the scaled chiral 
condensate is always a monotonic decreasing function of m. These points also qualitatively 
coincide with the lattice results [36, 37]. For T = 180MeV, pion superfluidity phase occurs at 
Hi = 75MeV which is very close to fi c and the phase transition is second order. The scaled 
pion condensate first increases with the increasing and then decreases with the increasing \xi 
until it abruptly vanishes at [i\ = 456MeV. With increasing T, the first order phase transition 
for the pion condensate vanishing becomes more and more soft and it evolves as a second order 
phase transition when T > 210MeV. Fig. 3 also shows that the two phase transition points 
for ^3 symmetry breaking and restoration get more and more close to each other with the 
increasing T and eventually coincide at ~(223MeV, 135MeV). Numerical results also suggest 
that the peak of dQ/djJLi is located at the vicinity of the pion condensate vanishing point. 

The reason that the pion condensate first increases and then decreases with the increasing \xi 
at fixed T is that it eventually becomes energetically favourable to produce fermionic excitations 

1 In this paper the Polyakov loop susceptibility is defined as the partial differential of 3> with respect to T or 
fii. The definitions of the susceptibilities for chiral condensate and pion condensate are alike. 
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which then competes with the production of fermions in bound states with anti-fermions, that 
is pions. It is reasonable that to produce quark excitations gets more and more advantaged 
in contrast to the production of quarks in pions with anti-quarks for high T, as suggested by 
Fig.3. Since this phenomena is also observed at T ~ 223MeV where the non vanishing pion 
condensate(with very small value) only appears in the vicinity of fxi ~ 135MeV, we believe 
that the above result is not an effect due to the ultraviolet cutoff. 

The (T,hj) phase diagram of two flavor PNJL model is shown in Fig. 4(a), in comparison 
with the standard NJL phase diagram described by Fig. 4(b). We only plots the phase diagrams 
in the range [i\ < 500MeV since the investigation for more high isospin chemical potential is 
beyond the capability of the NJL type models. The PNJL phase diagram is qualitatively 
consistent with the lattice results given in [37]. Fig. 4(a) shows that the pion condensate 
always vanishes above T c °. The TCP is located at (210MeV, 300MeV), which means that the 
phase transition for I3 symmetry restoration in the region \ii > 300MeV is always first order. In 
the rough, for fij > 135MeV, the chiral condensate is greatly suppressed and the melting of pion 
condensate is simultaneously accompanied by the crossover phase transition for deconfinement; 
while for \i c < m < 135MeV, the chiral condensate still keeps a considerable value when the 
pion condensate vanishes and the crossover phase transitions for chiral symmetry restoration 
and deconfinement are still coincide. Note that numerical study suggests that the standard 
NJL model also sustain the existence of the TCP in the (T,fij) phase diagram. The TCP 
obtained in standard NJL is at the vicinity of (140MeV, 400MeV). However, the standard NJL 
has no capability to give the information concerning the deconfinement phase transition. 

3.2 The case for A ^ and fi B = 

The previous lattice simulations for pion superfluidity phase transition are performed by intro- 
ducing an explicit isospin symmetry term which is proportional to a small parameter A in the 
lagrangian. The phase diagram is then obtained by extrapolating the results from A ^ to 
A = 0. To compare with the lattice data, we also investigate the A / within PNJL formalism. 
In addition, studying A / case can further deepen our understanding the role of the Polyakov 
loop dynamics . 

with A / 0, the pion condensate is no longer the true order parameter and the phase transition 
for pion superfluidity becomes a crossover. We perform the calculations with A = 5MeV(with 
the same order as mo), 0.5MeV and 0.05MeV, respectively. Fig. 5 gives the scaled pion 
condensate, chiral condensate and Polyakov loop as functions of isospin chemical potential at 
[is = and T = 210MeV for different A. Fig. 5(a) shows that when A = mo the chiral 
condensate and pion condensate will have the same magnitude in vacuum, which indicates that 
a small explicit symmetry term has significant influence on the vacuum structure. With A 
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decreasing, the magnitude of n in m < fi c is greatly suppressed while its influence on the large 
Hi region is not so significantly. We can see that Fig. 5(c) which is obtained with A = 0.05MeV 
is almost identical with Fig. 3(b) which is obtained with A = 0. Fig. 5 also shows varying of 
small A has little impact on the shape of Polyakov loop <3?. Here we can get conclusion that the 
method taken in lattice simulation by extrapolating the results from A / to A = is reliable. 

For small A, we can use the peak of the pion condensate susceptibility to indicate the crossover 
phase transition for /3 symmetry breaking or restoration. In Fig. 6, we present dir/d^ii and 
d^/dfii obtained with the same parameters as Fig. 5. The curve of dir/d^i has two peaks, 
which correspond to the crossover phase transitions for .Z3 symmetry breaking (the left peak) 
and ^3 symmetry restoration (the right peak), respectively. It is shown that the peak of d^/dfii 
perfectly coincides with the right peak of dir/dfii. Fig. 6 also shows that the peak position of 
d^/dfii is insensitive to the value of A when A << tjiq while the extent of superposition between 
the peaks of d^/dfii and dir/dpi is closely related to the magnitude of the explicit symmetry 
breaking term. In comparison with the Fig. 3(b), we find that the peak position of crossover 
phase transition for 1% symmetry breaking or restoration is considerably consistent with the 
corresponding true phase transition point for ^3 symmetry breaking or restoration obtained 
with A = 0. This results confirm that the melting of pion condensate for > fi' c « 135MeV is 
simultaneously accompanied by the crossover deconfinement phase transition, just as shown in 
Fig. 4(a). 

The evidence for the existence of TCP is illustrated in Fig. 7. Fixing A = 0.5MeV << rriQ, 
for high enough temperature(T > 210MeV), the crossover phase transition for ^3 symmetry 
breaking is similar to the crossover phase transition for ^3 symmetry restoration, which in- 
dicates that both phase transitions are second orders for A = 0; While for the low enough 
temperature(T < 210MeV), the right peak becomes significantly steeper than the left one, 
which indicates that it corresponds to a first order phase transition for A = 0. Fig. 7 also 
shows that the two peaks get more and more close to each other with increasing T and eventu- 
ally coincide at a point (T', //). For A << mo, we find that the point (T 1 ',//) is almost identical 
with the point (223MeV, 135MeV). 

Above results can also be confirmed by exploring the peak positions of dn/dT and <9<3?/<9T, 
which is shown in Fig. 8. For low fij = 80MeV, we can see that differences among three 
crossover phase transitions for chiral condensate, pion condensate and Polyakov loop are still 
within a few MeV with A « uiq. For the same A, the peaks related to the pion condensate and 
Polyakov loop both get more and more steep with increasing m, which implies that second 
order phase transition for ^3 symmetry restoration eventually involves into the first order phase 
transition. 
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3.3 Comparison with Lattice results 

For quantitative comparison with the corresponding lattice data, following [28], we also reduce 
the Tf} by rescaling the parameter To from 270MeV to 190MeV. In this case, though the 
perfect coincidence of the crossovers for the deconfinement and chiral phase transitions is lost, 
the difference between the two transition is within 20MeV[28]. The average of these two 
temperatures can be defined as T c ° with the value 180MeV, which is in the range of recent 
lattice results [46, 47]. Note that in the following calculations, we keep the liquid quark mass 
varying with the change of (3 and other parameters of the PNJL unchanged. 

The PNJL results are compared with the lattice data obtained at fixed m = 0.05, fij = 0.4 
and A = 0.005 with N t = 4[36] in Fig.9(the PNJL results were calculated with the same 
parameters). The critical (3 C at zero chemical potential is 5.3198(9) for m=0.05 [37], which may 
corresponds to T® = 180MeV. In lattice simulations with fixed Nt, the temperature T varies 
with lattice spacing a and T = N t a~ l . Using 2-loop running of the coupling g(a), the lattice 
spacing a can be expressed as 

-AZ'l^f^expf-JL), (22) 

where Q = 6/j*. A = ^(Uf - ??) and ft = ^(2^ - « - &^-). The 
one dimensional constant parameter Al can be determined by comparing the known critical 
f3 c with the corresponding critical temperature T®. Using this relation, we can transform Ts 
into f3s and compare the PNJL results with the corresponding lattice data. Note that in lattice 
simulations, every quantity is made dimensionless by multiplying the appropriate power of 
a(J3). 

Fig.9(a) shows that the dependence of Polyakov Loop on [3 in PNJL is well consistent with 
the lattice data. The peak of Polyakov Loop susceptibility from PNJL is located at (3 = 5.27, 
which coincides with the lattice result depicted in Fig. 14 in Ref.[36]. With 2-loop running of the 
coupling g, this crossover transition point for deconfinement corresponds to T = 169MeV and 
Hi = 270MeV(note that m also changes with (3) with m = 34MeV. The variation of the pion 
condensate with respect to (3 is shown in Fig. 9(b). Since the original data for pion condensate 
in [36] are all dimensionless quantities scaled by a series of physical units (corresponding to 
different lattice spacing a{(3)), we presented these data in physical units and then rescaled the 
obtained results by the pion condensate at (3 = 5.2. Both PNJL and lattice data show that the 
pion condensate first increases with (3 and then decreases with increasing (3. The crossover phase 
transition for I3 symmetry restoration in PNJL model takes place at (3 = 5.3 or (T = 175MeV, 
III = 280MeV) with m = 35MeV. In contrast with the critical value shown in Fig.9(a), both 
the temperature and isospin chemical potential of this point are shifted by less than lOMeV. 
We also performed calculations in PNJL for m = 0.05 by keeping Tq = 270MeV and adopting 
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T c = 222MeV as in previous sections and found that both crossover transition points located at 
the same point (3 = 5.3. Therefore the coincidence of the crossover transitions for deconfmement 
and ^3 symmetry restoration obtained in lattice simulation [36] can be successfully interpreted 
by including the Polyakov Loop dynamics in the NJL model. 

The 16 3 x 4 lattice simulations in [37] with the set of same parameters as in Fig.9 strongly 
suggest that there would be a first order phase transition for the ^3 symmetry restoration 
between (3 = 0.525 and (3 = 5.3 in the A — ► limit. The authors of [37] claimed that this point 
is very close to the TCP. This first order phase transition was confirmed in our calculations 
within PNJL for A = and the transition point is located at f3 = 5.32 or (T = 180MeV, 
Hi = 288MeV), which is also in the vicinity of the TCP obtained in PNJL (m = 36MeV). 

At fixed f3 = 5.0 or T = 119MeV, the pion condensate as a function of fij for m = 0.05 or 
m= 24MeV and different A is plotted in Fig. 10. The range of m in this plot corresponds to 
< m < 239MeV. The ellipses are the linear extrapolation of the lattice data to A = and 
the solid line corresponds to the PNJL result obtained at A = 0.0005. Fig.10 indicates that 
the PNJL results are quite well consistent with lattice data for A = 0.01, 0.005. This plot also 
shows that the linear extrapolation to A = in lattice simulation is fairly good agreement with 
the PNJL data in the low and high \xj region, while the deviations are explicit in the vicinity of 
the true phase transition point. In contrast to the critical point m = 0.25(119MeV) obtained 
from the lattice linear extrapolation, the critical isospin chemical potential in PNJL model 
is m = 0.3(143MeV). Both results sustain that the I3 symmetry broken phase transition is 
second order for low temperature (in contrast to T®) and no decrease of pion condensate with 
increase of fij is observed in the plotted range of fij . 

4 Conclusion and outlook 

The PNJL model represents a minimal synthesis of the two basic principles which govern 
QCD at low temperatures: spontaneous chiral symmetry breaking and confinement. In this 
work, we extended the two flavor Polyokov-loop-extentded NJL model to the finite isospin 
chemical potential case, which can be investigated by lattice simulations. By including the 
pion condensate degrees of freedom, the trace of Polyakov loop and its conjugate can still 
be extracted out and the formalism which can simultaneously couple the fields of the pion 
condensate, chiral condensate and Polyakov loop is obtained. 

Within this framework, the dependence of the pion condensate, chiral condensate and Polyakov 
loop on temperature and isospin chemical potential are explored and the two flavor (T, phase 
diagram is plotted at the mean field level. Our calculations confirmed that there exists a tri- 
critical point in the (T, fj,i) phase diagram. In the low fij region with pion condensate, the 
crossover phase transitions for chiral symmetry and deconfmement are still perfectly coincide 
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and the isospin symmetry restoration phase transition for high isospin chemical potential region 
is simultaneously accompanied with the crossover phase transition for deconfinement. We also 
directly compared the PNJL results with the corresponding lattice data at finite temperatures. 
In general, a qualitative agreement between these two methods is confirmed and in some cases 
even a quantitative agreement is obtained . Therefore, our results provides further test on the 
validity of this enhanced quasi-particle model. At the same time, the conclusions obtained 
within PNJL model also give a way to test the reliability of the lattice data for finite isospin 
density. 

In this paper, we only explored the (T, phase structure of two flavor PNJL at zero 
baryon chemical potential. In our forthcoming paper, we will investigate the Polyakov loop 
dynamics on the phase structure, thermodynamic and meson properties at finite temperature, 
isospin chemical potential and baryon chemical potential. Further exploration on the phase 
diagram and meson properties of three flavor PNJL model by including the strangeness chemical 
potential is also underway. 
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Figure 1: Scaled pion condensate, chiral condensate and Polyakov loop <5 as functions of 
temperature at zero baryon chemical potential for different fj,i with A = 0. 
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Figure 9: Polyakov Loop and pion condensate as a function of (3. The lattice data (the 
black circles) are taken from [36] for m = 0.05, A = 0.005 on an 8 3 x 4 lattice. In the right 
panel, the white circles correspond to the original lattice data which were first represented 
in the same physical units and then rescaled by the obtained pion condensate at (3 = 
5.2(vr'). 
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SU(3) N f =2 /3=5.0 m=0.05 




Figure 10: Scaled pion condensate as function of fij for A = 0.01, A = 0.005 and A— > 
at fixed m=0.05 and (3 = 5.0. The lattice data are taken from [36] and we rescaled them 
by the lattice result for the pion condensate at 2\xj = 1 and A = 0.01. The curves for 
the PNJL data are also scaled by the pion condensate obtained in PNJL at 2/ij = 1 and 
A = 0.01. 
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